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ABSTRACT 


Microgrid energy systems are one of suitable solutions to the available 
problems in power systems such as energy losses, and resiliency issues. 
Local generation by these energy systems can reduce the role of the upstream 
network, which is a challenge in risky conditions. Also, uncertain behavior of 
electricity consumers and generating units can make the optimization 
problems sophisticated. So, uncertainty modeling seems to be necessary. 
In this paper, in order to model the uncertainty of generation of photovoltaic 
systems, a scenario-based model is used, while the robust optimization 
method is used to study the uncertainty of load. Moreover, the stochastic 
scheduling is performed to model the uncertain nature of renewable 
generation units. Time-of—use rates of demand response program (DRP) is 
also utilized to improve the system economic performance in different 
operating conditions. Studied problem is modeled using a mixed-integer 


linear programming (MILP). The general algebraic modeling system 
(GAMS) package is used to solve the proposed problem. A sample microgrid 
is studied and the results with DRP and without DRP are compared. It is 
shown that same robustness is achieved with a lower increase in the 
operation cost using DRP. 
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1, INTRODUCTION 

Microgrid is a scaled-down power system that can benefit from local generation units to supply the 
energy demands. Optimal operation of microgrid is challenged by various factors including uncertainties, 
which can threaten the normal operation of microgrid energy systems [1-5]. Therefore, the robust operation 
of microgrid should be investigated. 

A brief summary of studied papers about microgrid opera-tions is presented in the following: 
Implementing a two-stage adaptive robust optimization-based collaborative operation, the economic and 
collabo-rative operation of multi-microgrids has been evaluated as a unit commitment problem [6]. The total 
operation cost of multi-microgrids was reduced subject to uncertainty of so-lar power systems. In [7], optimal 
design of a hybrid electric power generation system for isolated zones has been inves-tigated and the total 
cost of power generation was mini-mized using particle swarm optimization. With the aim of calculating the 
sharing power of distributed energy re-sources with uncertainty and nonlinear load, the robust multi-objective 
control method has been employed in [8]. Optimal design and scheduling of power dispatch in a mi-crogrid 
has been achieved using a sliding time window op-timization modeling considering the unbalance of 
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biomass-integrated renewable energy and energy storage [9]. In or-der to optimize the performance of energy 
management in the presence of renewable energy resources, a multi-objective generation scheduling model 
was developed un-der demand response program [10]. By implementing de-mand response program (DRP), 
the optimum components sizes in a microgrid was evaluated to decrease the cost of energy in [11]. With the 
aim of solving microgrids operation problem in the presence of electric vehicles, the robust op-timization has 
been implemented to reduce the operation cost of microgrid in [12]. In order to improve the characteris-tics 
of electric power distribution systems, conventional power generators and various types of renewable energy 
resources and electric vehicles have been used to improve the voltage profile and decrease the cost of power 
losses and energy [13]. Robust technique has been employed to control an isolated microgrid in the presence 
of high pene-tration of renewable energy sources [14]. Risk-based optimal operation of microgrid system is 
evaluated using robust optimization method in the presence of electric vehicles in [15]. Distribution system 
with multi-microgrid is optimally scheduled under uncertainty with robust optimization method in [16]. 
Similarly, robust optimization approach 1s used in [17] to solve uncertainty-based optimal energy transactions 
problem in multi-microgrid systems. In [18], a brief summary of studied papers about operation of pho- 
tovoltaic based microgrids has been presented to investi-gate various combination methods of energy storage 
and electric vehicles with other energy resources. It is notewor-thy that optimal scheduling of microgrid 
systems have been studied under different types of demand response programs in several research papers 
[19-21] which repre-sents essential role of these programs in improvement of optimization results. 

In this paper, the robust operation of microgrid energy system subject to uncertainties of load and 
renewable gen-eration unit is studied considering the time-of-use rates of demand response program. In order 
to model the uncer-tainty of generation of photovoltaic systems, a scenario-based model is used, while the 
robust optimization method is used to study the uncertainty of load. Combination of robust optimization 
method, which is a powerful strategy providing uncertain modeling method, with stochastic programming to 
model uncertainty ensures optimal reaction of operator against possible negative consequences of uncer- 
tainty. Furthermore, demand response can enhance effec-tiveness of these methods in uncertain handling. 

It should be noted that the level of robustness depends on the policies taken by the operator of 
microgrids. According to each desired level of robustness, operation cost of mi-crogrid is increased. 
Therefore, robust optimization method can be used as a powerful uncertainty modelling tool to provide 
operation strategies for energy systems against different uncertainties. The studied problem is schematical-ly 
illustrated in Figure 1. 
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Figure 1. Studied problem 


1.1. Paper Structure 

The remaining parts of the paper are organized as: Robust operation of microgrid energy system is 
formulated in Sec-tion 2. The case studies and results of the proposed robust optimization scheme are 
presented in Section 3. Finally, the paper is concluded in Section 4. 


2. PROBLEM FORMULATION 
In this section, robust operation of microgrid energy system with load uncertainty is modelled, 
considering the demand response program. 


2.1. Objective Function 


As the objective function, the cost of the purchased power from the upstream power grid (1) should 
be minimized. 


T 
Min Cost = >, A, x PR" 
t=1 (1) 
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net 


In (1), Cost is the total operation cost of the microgrid; A, and f are the price and the amount of 
the power purchased from the upstream network. 


2.2. Constrains 
The power balance is formulated as 


Load, DRP Ess,ch PV t Ess ,di 
| tai ae a =f sl Shek 


(2) 
Load ,DRP p®s ch 
where “' is the load with DRP; ~' is the charging power of the battery energy storage systems 
PV net 
(ESS); e is the power generated by the solar system; ~‘ is the power from the upstream network; and 
Ess,dis 
Z is the discharging power of battery ESS. 
The active power generated by the photovoltaic system is expressed as 
PE = AXR,XMpy (3) 
where A is the installation area: K, is the solar radiation; and /?v is the PV output efficiency. 
The power imported from the upstream network is limited as 
ates age ote 4) 


net net 


Pp. ie 
where ™* and ™" are the maximum and minimum allowable powers. 
The stored energy in battery ESS depends on the state of charge of the battery at previous time step 
and the charging/discharging power at current time interval. The available energy is 


Ess __ Ess Ess,ch ch Ess,dis dis 
Wawa Pe gy pet fy 


(5) 
W Ess Ess,ch Ess,dis 7 ch 
where ‘is the stored energy level; “‘ and “' are the charging and discharging powers; and 
dis 
and ” are the charging and discharging efficiencies of the battery ESS, respectively. 
The stored energy of the battery ESS is limited as 
Ess Ess Ess 
Wain We S Wrrax (6) 
We wes 
where ™“ and ™” are maximum and minimum energy of battery ESS, respectively. 
The charging and discharging powers of battery ESS are limited by (7)-(8). 
a“ epee < Boe < oF ne eae (7) 
a xP < Bee < a ny cel (8) 
pee pas pee pee 
Where  ™® max min’ and © ™* ~—s are the minimum and maximum discharge and charging power 


limits of the battery ESS, respectively. 
j ch 


dis 
In (7) and (8), “and @ are binary variables of discharging and charging of the battery ESS, 
which are limited as (9) since the battery ESS 1s either in charging or discharging. 


ch dis 
a +a <1 
; (9) 


The demand response program with the time-of-use method is modeled based on (10)-(13). 


‘pice ,DRP = poe 


shup shdo 
t t 


ep =D (10) 
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Load shup shdo 
where “‘ 1s the base demand, Pi and + — are the increased and decreased values of the base load in 


DRP, which are limited as (11) and (12), respectively. 


0 < Do < LPF’ x poe x | bad 


(11) 
shdo shdo Load shdo 
O< p< LPF’™ xP™ xI (12) 
shup shdo a ; : I shup ] shdo 
where LPF and LPF’ are the limitation of the increase and decrease in load, and ‘ and‘ are 
binary variables to represent the increase of load in DRP. 
The increased and decreased powers for DRP are balanced as: 
- h S hd. 
pa a a a 
(13) 


where / is the scheduling horizon. 

The robust optimization method and the different uncertainties are modelled as (14)-(20) to 
formulate the robust operation of microgrid energy system [22-23]. As shown, after implementing robust 
optimization, deviation of uncertain parameter are added to the objective function (14). (15) refers to the 
constants of base problem. As shown in (16)-(20) are used to model deviation of uncertain parameter. 
Then, by applying formulations presented in [22-23], the studied problem under uncertainty can be 
expressed as: 


Min (S A, x PR + 20) + 3 tu 


a a (14) 
S.C. 
(2)-C13) (15) 
Zy+q,,24,y., t=1,....T 
(16) 
G20. Faihegt (17) 
y, 20, t=1,...,7 (18) 
Z 20 (19) 
eee ia (20) 


It should be noted that “° and 4 are dual variables of optimization problem. Duel variables are 


formed when dual form of an optimization problem is written using duality theorem. Also, Yo is an integer 


defined to determine either uncertainty is considered or not. When uncertainty is considered, 'o=! 
T, =0 


and 
otherwise 


3. NUMERICAL STUDIES 

Robust operation of microgrid energy system is evaluated considering the uncertainties and the DRP 
in this section. The input data used for simulations is presented in Subsection 3.1. The proposed robust 
optimization problem is solved and the results without and with DRP are given in Subsections 3.2 and 3.3 in 
the following, respectively. The studied test system is captured in Figure 2. 
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Figure 2. Test systems 


3.1. Input Data 


The hourly solar radiation is illustrated in Figure 3 [24]. The price of the power imported from the 
upstream network is illustrated in Figure 4 [25]. The energy demand of microgrid in different levels is 
illustrated in Figure 5 [25]. As illustrated, load is assumed to fluctuate between 0.7 and 1.3 of its base value. 
The technical data of battery ESS 1s provided in Table |. It should be noted large number of scenarios are 


generated and then reduced for renewable units and the expected values are calculated to be used in the 
optimiza-tion process. 


Table 1. ESS Characteristics 
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Figure 5. Upstream network price 


3.2. Results of Robust Energy Management without DRP 


In this section, results of the robust and stochastic scheduling of the microgrid energy system are 
presented. Calculating the expected values of solar radiation in each hour, the robust optimization method is 
applied to provide strategy for different levels of load uncertainty without DRP. The load is assumed to 
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fluctuate between 0.7 and 1.3 of its base value. It should be noted that these values are selected random and 
can change based on the expectations of operator. Therefore, the most possible value of the load which is 
threatening the usual operation of the microgrid is 1.3 times the base value. Accordingly, the simulations are 
performed for several steps and the results are obtained. The total operation costs of microgrid against 
different levels of load uncertainty are presented in Table 2. 


Table 2. Robust Cost without DRP 
Load uncertainty level (of base value) 0.7 0.8 0.9 1.0 al 12 1.3 
Operation cost of microgrid (€) 829.779 956.269 1082.759 1209.249 1335.739 1462.229 1588.719 


According to the results presented in Table 2, in order to become robust against 30% variation in 
load, the operation cost of microgrid energy system should be increased by 31.38, to guarantee its normal 
operation against 30% additional load. On the other hand, when the load value is reduced to 70% of its base 
value, the operator of microgrid energy system can gain an economic profit up to 379.47 €. In fact, due to 
reduction of load value up to 30 %, daily operation cost of microgrid energy system will decrease up to 31.38 
%. Therefore, by employing each one of the provided strategies by robust optimization method, the energy 
system can be resilient against different levels of load uncertainties. As a result, the imported power from the 
upstream network in different load uncertainty levels without DRP is obtained which 1s illustrated 
in Figure 6. 
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Figure 6. Upstream network power without DRP 


In detail, when the uncertainty of load is not considered, the total power provided by the upstream 
network is 1881.753 kW. When the load is in maximum level, total imported power from the upstream 
network is 2466.453 kW. On the other hand, when the load is minimum, the total imported power is 
1297.053 kW. 


3.3. Results of Robust Energy Management with DRP 

Applying the time-of-use program of the DRP, the total operation costs against different levels of 
load uncertainty are obtained and presented in Table 3. It can be seen that the maximum value of load 
uncertainty can be handled by 5.05 % less increase in operation cost under DRP which is economically 
beneficial for the microgrid operator. Moreover, by reduction of load value to 70% of its base value, 
the microgrid energy system can benefit 9.68% more reduction in cost utilizing the DRP. The power 
imported from the upstream network in this case is illustrated in Figure 7. As shown, the load is shifted from 
the peak periods to off-peak ones, and thus the power is purchased economically optimal to reduce the total 
operation cost as much as possible. 


Table 3. Robust Cost with DRP 
Load uncertainty level (of base value) 0.7 0.8 0.9 1.0 ia) 12 1.3 
Operation cost of microgrid (€) 749.399 875.889 1002.379 1128.869 1255.359 1381.849 1508.339 
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Figure 7. Upstream network power with robust scheduling and DRP 


3.4. Comparisons 
To approve the positive role of DRP on the uncertainty-based performance of the microgrid, 
simulation results are presented in Table 4, with and without taking the DRP into consideration. 


Table 4. Comparison Results 


Load uncertainty level Operation cost of microgrid (€) 
(fraction of base value) Without DRP With DRP 
0.7 829.779 749.399 
0.8 956.269 875.889 
0.9 1082.759 1002.379 
1 1209.249 1128.869 
1.1 1335.739 1255.359 
1.2 1462.229 1381.849 
1.3 1588.719 1508 .339 


4. CONCLUSION 

In this paper, robust and stochastic operation of microgrid energy system is studied considering the 
uncertainties in load and renewable power of PV system. A large number of scenarios are generated first and 
then reduced using the fast-forward scenario reduction algorithm. Consequently, the expected values of 
uncertain parameters related to renewable power of PV system are calculated. The robust scheduling of the 
microgrid energy management is formulated as a mixed-integer linear programming (MILP) optimization 
problem. According to the results, robust scheduling of microgrid energy system in the worst condition is 
obtained against 30% higher load level compared to the base load, while daily operation cost of microgrid is 
increased by 31.38% without DRP. However, the same robustness can be achieved with a lower increase in 
the operation cost using the DRP, which is more economically beneficial for the microgrid operator. 
The proposed robust optimization method can be used as a powerful uncertainty modelling tool to provide 
operation strategies for energy systems against different uncertainties. 
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